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The mass spectra of ground-state triply heavy baryons are systematicaUy unscrambled and com- 
puted in QCD sum rules. With a tentative (QQ) — (Q') configuration for QQQ' , the interpolating 
currents representing the triply heavy baryons are proposed. Technically, contributions of the oper- 
ators up to dimension six are included in operator product expansion (OPE). The numerical results 
are presented in comparison with other theoretical predictions. 
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^ . I. INTRODUCTION 

1^ ' The triply heavy baryon, wherever light quarks are absent, is well and truly not a new topic but with 
a history. As one refers to the studies on their properties, it can be traced back to two more decades 
■ ago [l[ 01- However, contrasted with the singly and doubly heavy baryons (such as Refs. 0, 1^), only 
infrequent attention has been paid to the triply heavy baryons. Whereas, the case for triply heavy baryon 
may be improved and several approaches have already appeared in recent years, such as effective field 
theory lattice QCD [gl, variational method [7], bag model quark models [1, 10, 11 1 etc., for 



^ ^ which is gradually becoming a exciting and remarkable theme nowadays. First, the field of heavy hadron 
1^ _ spectroscopy is experiencing a rapid advancement, which is mainly propelled by the continuous discovery 
, ' of hadronic resonances (for reviews, e.g., see 13, l3|)- While experimentally reconstructing a candidate 



for ficcc is very difficult, it is not unthinkable according to Ref. 0. Especially for the startup of Largi 
Hadron Collider, it seems rather promising to establish triply heavy baryons in future lla l. Second, 

1^ — investigation of the triply heavy baryon is of great interest in understanding the dynamics of QCD at the 
0^ hadronic scale. Although the statement that QCD is the correct theory underlying strong interaction has 
. been commonly accepted and QCD is simple and elegant in its formulation, many questions concerning 
' dynamics of the quarks and gluons at large distances remain unanswered or, at most, understood only 
at a qualitative level. The quantitative description of the hadronic properties runs into however arduous 
0^ , difficulties. For example, it is a great challenge to extract information on the spectrum from the rather 
^P. ' simple Lagrangian of QCD. That's because low energy QCD involves a regime where it is futile to attempt 
. ^ perturbative calculations and, inevitably, one has to treat a genuinely strong field in nonperturbative 
1 methods. Briefly recapitulating the second point, triply heavy baryons, free of light quark contamination, 
are ideal prototypes to refine one's present understanding of heavy quark dynamics and may serve as a clean 
probe to the interplay between perturbative and nonperturbative QCD. Also stimulated by the above two 
aspects, it is interesting and significative to study their properties like masses through nonperturbative 
approaches, and the practitioner may resort to a vigorous and reliable working tool in hadron physics. 



the QCD sum rule 17|, which is a nonperturbative analytic formalism firmly entrenched in QCD. On 
the sum rule analysis, the triply heavy baryon systems are analogous to the cases of charmonium and 
bottomonium, where light quarks are all absent. In fact, there have alre ady been some works on calculating 
the charmonium and bottomonium masses in QCD sum rules, such as [18|. The c and b quark masses can 
also been determined from considering the two-point correlation function of the Q^^Q current (Q = c or 6), 
for instance in [TtI . [i9| , and some impressive progresses were made in updating the values of rric and m?, 
later, including the O(a^) perturbative corrections [i^. In addition, the semileptonic decays of have 



been investigated by three-point sum rules [2l|. Thereby, it is feasible for QCD sum rules to study triply 
heavy baryons and we would like to carry out the sum rule calculations of their spectra in this work. 
The paper is organized as follows. In Sec. |TT1 QCD sum rules for the triply heavy baryons are introduced. 
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and both the phenomenological representation and QCD side are derived, fohowed by the numerical analysis 
to extract the spectra and a comparison with other theoretical calculations in Sec. IIIII Section HVl contains 
a brief summary and outlook. 



II. TRIPLY HEAVY BARYON QCD SUM RULES 

A generic QCD sum rule calculation consists of three main ingredients: an approximate description of 
the correlator in terms of intermediate states through the dispersion relation, a description of the same 
correlation function in terms of QCD degrees of freedom via an OPE, and a procedure for matching these 
two descriptions and extracting the parameters that characterize the hadronic state of interest. Concretely, 
coming down to the mass sum rules for triply heavy baryon QQQ' (here Q and Q' can be the same or 
differently heavy quarks, c or 6), the starting point is the two-point correlation function 



n(g2)^j j d4xe'^-(0|r[j(x)j(0)]|0). (1) 

Lorentz covariance implies that the correlation function ([T]) has the form 

Ii{q^)^iIi^{q^)+Ii2{q^). (2) 

For each invariant function Hi and 112, a sum rule can be obtained. 

In the phenomenological side, the correlator can be expressed as a dispersion integral over a physical 
spectral function 



n(q2) = A|f ^ + - / '^''^ 2^ + subtractions, (3) 



where Mh denotes the mass of the triply heavy baryon. In obtaining the above expression, the Dirac and 
Rarita-Schwingcr spinor sum relations, 

Y,N{q,s)N{q,s) = i + MH, (4) 

S 

for spin-^ baryon, and 



H 



for spin- 1 baryon, have been used. 

In the OPE side, the correlation function can be written in terms of a dispersion relation as 

n.(g2)= r ds^^, z = l,2, (6) 



/(2mQ+mQ,)^ ^ 

where the spectral density is given by the imaginary part of the correlation function 

p,{s) = -Imn,?PE(,). (7) 

TT 

In detail, the spectral densities are calculated and embodied in Sec. Ill Bl 

After equating the two sides, assuming quark-hadron duality, and making a Borel transform, the sum 
rules can be written as 

Xle-'"-/^''' = r dsp^[s)e-^'^\ (8) 

J(2mQ+niQ/)2 
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X^Mne-''-/''^ = / dsp2{s)e'^/''\ (9) 

To eliminate the baryon coupling constant Xh , one reckons the ratio of derivative of the sum rule and itself 
and yields 

/•so , 2 ,2 

M% = I dspi{s)se-'/^^ I / dspi(s)e-^/^^ , (10) 

Ml = r dsp2{s)se-'/^'" / r dsp2{s)e~'/^'\ (11) 

A. The interpolating currents 

In a tentative picture for QQQ' system, the Q' orbits the bound QQ pair. The {QQ) — {Q') structure 
may be described similar to QQ' mesons, where the QQ pair plays the same role of the antiquark Q in 
QQ'. The study of such configuration can help one to adopt the appropriate interpolating currents. For 
the ground states, the currents are correlated with the spin-parity quantum numbers 0+ and 1+ for the 
heavy QQ diquark system, along with the quark Q' forming the state with = and the pair of 
degenerate states. For the latter case, the QQ diquark has spin 1, and the spin of the third quark is either 
parallel, = 1 , or antiparallel, = ^ , to the diquark. The choice of Ffc and F^, matrices in baryonic 
currents may be determined according to the rules in |22| . For the baryon with — i the current may 



be gained using SU{3) symmetry relations [23|. Consequently, following forms of currents are adopted 



JQqqq = £abc{QaCTkQb)^kQc, 

jn^^^, = eabc{QlCrkQb)r',Q',, (12) 
nQQ' = £abc^[2{QlCruQ',)T',,Q,. + {QlCTkQb)r'M 

^QQ' ^ eabc{Qlcr,Qt)r[Q',, 

Here the index T means matrix transposition, C is the charge conjugation matrix, a, b, and c are color 
indices, with Q and Q' denote heavy quarks. The categories of ground-state triply heavy baryons and the 
choice of F^ and F^, matrices are listed in TABLE HI 



TABLE I: The choice of F^ and Fj. matrices in baryonic currents. The index d in Sd, Ld, and J^"^ means diquark. 
{QQ} denotes the diquark in the axial vector state and [QQ] denotes diquark in the scalar state. 



Baryon 


quark content 




Sd 


Ld 


■^d 


Ffc 


Ffc 




{QQ}Q 


3 + 
2 


1 





1+ 


7m 


1 


^QQQ' 


{QQ}Q' 


1 + 
2 


1 





1+ 


7m 


7m 75 




{QQ}Q' 


3 + 
2 


1 





1 + 


7m 


1 


^'qqq' 


[QQIQ' 


1 + 
2 








0+ 


75 


1 



B. The spectra densities 

Implementing the calculation of the OPE side, we work at leading order in as and consider condensates 
up to dimension six. To keep the heavy-quark mass finite, one uses the momentum-space expression 
for the heavy-quark propagator. The final result is dimensionally regularized at -D = 4. It should be 



4 



distinguished for spectral densities of two sort triply heavy baryons, namely, containing the same heavy 
quark or differently. First, with 

o ramax 1-01 — nfi)'^ 

Pi{s) = da{J^ d(3 + d(3} y^^j; [ap{l-ap)s-{a^ + (3^ + ap + l)ml]s 



+ 



[a + /3)4 
a/3(l - a/3) 



j ""^ '^'^ij^' dp + j° d/3}^^^^[a/3(l - a/3)s - (a^ + /32 + a/3 + l)m^] 

a2/32(l _ a/3) 



da{ / d(]+ dl3} 



02 



[a + /3)2 



P2{s) = j"^^'' da{l\p + dp}^^^-^^[ap{l-ap)s-{a^+f3' + ap+l)ml]s 

+ ^ruQ J da{ d(3 + d/3}^-^[a/3(l - a/3)s - (a^ + 0^ + ap + l)m^]2 
+ / """^ '^"i '^Z? + /° '^/5}r^^[«/3(l - - + /3' + a/3 + l)m|] 



for Oqqq. The integration limits are given by 



(a + /3)^ 



\ 



(s2 - 6m|s - 3m^) - (s - m^)^ 


/is-ml){s 


- ^^D 




(s2 - 6to^s - 3m^) + (s - m|)^ 


/is-ml){s 


-9m^) 



a{s — mg) — JoP's'^ — Gmga^g „ 4771^ _ s^^a^ — 4mga"^s 
/3i = ;Tr^3 — ^ T\ ^and 



2(a2s + 77i|) 



/32 



Next, with 



i(s — TOg) + ^J a^s^ — GTTiga^s — 4mg — 3TOgQ;2 — irriQa^i 
2(a2s + m|) 

a2/32(l-a/3)2 



o ramax rPi /■- rv^/^-^fl — n//3r 

= £ '^"i / + ^/^} (a + /3)4 ["^(^ - - + + (1 - "/^)"^Q')]« 

52 /-ct 



2^4 

3 



/3i 



da{ I d(3+ dp} 

J 02 

r0 

dp + 

10 



{a + /3)4 
a/3(l - a/3) 



[a/3(l - ap)s - {{a + pfm^ + (1 - a/3)m^0]' 



(a + /3)4 

/"max /-Pi /"t q;/3 

da{ dp + dP}j^-^^[ap{l - ap)s - ((a + pfml + (1 - a/3)m|0 

257r4 



/ da{ dP+ dP}^ , 

da{/ 'd/3+ ["dP} 

Jo J 02 + P) 



CXmin 
^max 



"^^^ "f' [a/?(l - a/3)s - ((a + pfml + (1 - a/3)m^,)]s 



27r4 



da{ r dp + 1^ dp}^—^[aP{l - aP)s - {{a + pfm^ + (1 - a/3)m|,)]' 
-'0 J 02 I" + Pi 



3 2 



IqITIQ' 



/ da{ df3+ d(3}-——[a0{l - af3)s - {{a + pfml + (1 - a/3)m^,)] 

Ja„,,i„, Jo J 02 [a-tp) 



{g^G^} r"'"" , ( f .o p (l-a/3)[(l-a/3)2-a/3(a + /3)2 



for Oqqq/ , 



+ {a^ + l)][a/3(l - ap)s - {{a + pfml + (1 - ap)ml,)]s 

+ ^ da{ dfi + r dPj .'f^-fl^M a + 

2^ ia^in Jo J 02 (a + /?f(a^ + l) 

+ (a^ + l)][ap{l - ap)s - {{a + pfml + (1 - a/3)m^0]' 

I ramax r0i ^ 

+ ^ X^^^ I + (^^TWTT) [4(1 - + 

+ I3{a^ + l)m^] [a/3(l - a(3)s - ((a + (3fml + (1 - a/3)m|,)] 
(52G2) „.2/ 



+ 



132 

+ (1 - a/3)(a2 + l)mQ,][a(3{l - a(3)s - ((a + (3fm% + (1 - a(3)m%,] 

1 ^ ^ 

+ (1 - a/?)(a2 + l)mQ-][a/3(l - a/3)s - ((a + /3)2m^ + (1 - a/3)m^,)]' 



1 (-amaoi rfii f-t I 



+ (a^ + l)][ap{l - ap)s - {{a + (if ml + (1 - a/3)m^, 

(ff^G^) ^ /'^^ , ^^^^ 1 

^ 3 • 2-^4 ^ 

+ (1 - al3){a^ + 1)[(1 - a/3)2 - + /3)2]mQ/}, 



■ da{ d/3 + d/3}----i__-{4/3(a + I3f[0'{a + /?) - a(l - ap)]mQ 

Jarr.,r, ^ J 02 (« + /3) ^ + 1 ) 



foi' f^QQQ'' and 

Pi{s) = -^J^ da{ I dp + d(3} (J^^^4 M(l - - ((« + + (1 - aP)ml,)]s 

- ^ da{f^^ df3 + d/3}^gi-^[a/3(l - a/3)s - ((a + (if ml + (1 - a/3)m^0]' 

- da{j^^' d(3 + ^° dP}j£^[afi{l - afi)s - ((a + (if ml + (1 - a/3)m^0] 

2^^ io -//Ss (Q! + /3)^ 

P2{s) = -^niQ, da{j^^^ dp + rf/3}^fc^[a/3(l - a/3)s - ((a + /3)^m^ + (1 - aP)rnl,)]s 

O rOimax r0i fk 1 _ ^« 

mQ, / da{ dp+ dp}-—-^[ap{l-ap)s-{{a + pfml + {l-ap)ml,)f 

Joimir^ Jo J 02 + 



2%4 
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247r' 
267r4 



da{/ d(3+ d(3}-—- 
da{ df3+ [" d(3} 

Jo 



a/3{l - a/3)s - {{a + /3fm^g + (1 - a/3)'TiQ/)] 
(1 - ap)[{l - + apja + 



for Q'qqq,. The integration limits are given by 



Pi 

P2 



\ 



(s2 - 4to^s - 2mLs - 4:m%m%, + ml,) - (s - "iL)^[s - (2mQ - toq')2][s - (2toq + mQ'^] 



SrriQS 



(s2 _ 4to^s - 2m^,s - 4TO^m^, + m^,) + (s - "^q') VI'S - (2toq - toq')2][s - (2mQ + mQ/)^] 



a{s — 2mQ + niQ,) — ^a'^[s — {^ttIq — TOq/)]^ — 4{a'^mQ + m'Q,){m'Q + a^s) 

2(a2s + m|) ' 

a(s - 2m| + m|,) + y'a2[s - (2m^ - m|,)]2 - 4(a2m| + m^,)(m^ + a2s) 

2(a2s + 77i|) 



and 



III. NUMERICAL ANALYSIS 



In numerical analysis, not both but only the sum rule (jlOp will be numerically analyzed for brevity, and 
the input values are taken as rric — 1.25 GeV,m.b — 4.20 GeV with (g^G^) — 0.88 GeV*. Complying with 
the standard procedure of sum rule method, the threshold sq and Borel parameter are varied to find 
the optimal stability window, in which the perturbative contribution should be larger than the condensate 
contributions while the pole contribution larger than continuum contribution. Accordingly, the sum rule 



66b, 



windows are taken as ^/so = 5.5 
M2 = 14.0 - 15.5 GeV^ for 17, 
9.5 GeV, AP = 8.5 
yi^^= 11.5 - 12.5 GeV, AP 



6.5 GeV, A/2 = 6.5 - 8.0 GeV^ for fl 
= 8.0 - 9.0 GeV, AP = 8.5 

'so 



14.0 



10.0 GeV^ for fl^ 



10.0 GeV^ for f^*^^, ^ = 11.0 - 12.0 GeV, AI^ = 10.0 ~ 11.5 GeV^ for n. 



15.0 GeV, 
so = 8.5 ~ 

bbci 



for ^ccb^ and 



10.0 - 11.5 GeV^ for 17*^,^, ^ = 8.5 - 9.5 GeV, AP = 8.5 - 10.0 GeV" 
^ = 11.5 - 12.5 GeV, AP = 10.0 - 11.5 GeV^ for n^^^, respectively. The corresponding 
Borel curves are shown in Figs. 1-4. In Table [TTl the numerical results are presented, together with the 
predictions from other theoretical approaches. It is worth noting that uncertainty in our results are merely 
owing to the sum rule windows, not involving the ones rooting in the variation of the quark masses and 
QCD parameters. Note that the QCD 0{as) corrections are not included in this work, whose calculations 
for triply heavy baryons are quite complicated and tedious as one has to tackle some three-loop massive 
propagator diagrams. However, it is expected that the QCD 0{as) corrections might be under control 
since a partial cancellation occurs in the ratio obtaining the mass sum rules ([TU| and ([Tl]) . This has been 
proved to be true in the analysis for the singly heavy baryons (the radiative corrections to the perturbative 
terms increase the calculated baryon masses by about 10%) [24| and for the heavy mesons (the value of 
Jd increases by 12% after the inclusion of the 0{as) correction) 25[. After a detailed comparison, one 
can find that the central values of our results are lower than other predictions from potential models, in 
particular, for Qbtb, slightly more than 1 GeV, whereas the relative discrepancy approximates to 10%. In 
addition, for flccc, our result is in good agreement with the lattice QCD simulation in Ref. [6], but the 
other comparisons can not be made for lack of lattice results at present. 
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5.5GeV 
-e.OGeV 
-6.5GeV 




14.0GeV 
-14.5GeV 
-15.0GeV 



1 7.5 8 8.5 



3.5 14 1 4.5 M'(GeV') 15 15.5 16 



FIG. 1: The dependence on M for the masses of Vlccc and Q.hhh from sum rule jTUJl. The continuum thresholds are 
taken as = 5.5 - 6.5 GeV and ^ = 14.0 - 15.0 GeV. 



a.OGeV 
-a.SGeV 
-g.OGeV 



8.5 9 M'(GeV') 9.5 10 10.5 



8.5GeV 
-g.OGeV 
-9.5GeV 



9.5 10 10.5 



FIG. 2: The dependence on for the masses of flccb and ^^ct from sum rule (|10|l . The continuum thresholds are 
taken as ^ = 8.0 - 9.0 GeV and ^ = 8.5 - 9.5 GeV. 



II.OGeV 
-II.SGeV 
-12.0GeV 



9.5 10 10.5M^(GeV')11 11,5 12 



II.SGeV 
-12.0GeV 
-12.5GeV 



9.5 10 10.5M=[GeV=)11 11.5 12 



FIG. 3: The dependence on A/ for the masses of flbbc and $1^^^ from sum rule (|10|l . The continuum thresholds are 
taken as ^ = 11.0 - 12.0 GeV and ^ = 11.5 - 12.5 GeV. 




8.5 9 M'(GeV') 9.5 10 10.5 



II.SGeV 
-12.0GeV 
-12.SGeV 



9.5 10 lO.SM^iGev'lll 11.5 12 



FIG. 4: The dependence on M for the masses of ^^cb ^^^^ ^bbc froni sum rule (|10p. The continuum thresholds are 
taken as ^ = 8.5 - 9.5 GeV and ^ = 11.5 ^ 12.5 GeV. 
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TABLE II: The mass spectra of triply heavy baryons (mass in unit of GeV). 



Baryon quark content 


Sd 


Ld 


'^d 


This work 


Ref. \1\ 


Ref. [2] 


Ref. [6] Ref. [7] Ref. [8] 


Ref [101 


" "CCC 


{cc}c 


3 + 
2 


1 





1+ 


4.67 ± 0.15 


4.79 


4.925 


A fifti A 
4.Doi 4. / D 


A 777 


4.803 


\ Ihhh 


{bb}b 


3 + 
2 


1 





1+ 


13.28 ± 0.10 


14.30 


14.760 


14.37 


14.276 


14.569 


^ccb 


{cc}b 


1 + 

2 


1 





1+ 


7.41 ± 0.13 








7.984 


8.018 


^ccb 


{cc}b 


3 + 
2 


1 





1+ 


7.45 ±0.16 


8.03 


8.200 


7.98 


8.005 


8.025 


^bbc 


{bb}c 


1 + 

2 


1 





1+ 


10.30 ±0.10 








11.139 


11.280 


^bbc 


{bb}c 


3 + 
2 


1 





1+ 


10.54 ±0.11 


11.20 


11.480 


11.19 


11.163 


11.287 


^ccb 


[cc]b 


1 + 
2 








0+ 


7.49 ±0.10 












^bbc 


[bb]c 


1 + 

2 








0+ 


10.35 ±0.07 













IV. SUMMARY AND OUTLOOK 

In a tentative {QQ) — (Q') configuration, the QCD sum rules have been employed to calculate the 
masses of triply heavy baryon QQQ' {Q = Q' oi Q Q'), including the contributions of the operators 
up to dimension six in OPE. The mass values extracted from the sum rules are collected in comparison 
with other theoretical predictions. The results in this work are lower than the predictions from potential 
models, nevertheless, the one for flccc is well compatible with the existing lattice study. Indubitably, there 
are still plenty of problems desiderated to resolve. Experimentally, the evidence on triply heavy baryons 
are expected to reveal nature of them, especially after the putting into operation of the Large Hadron 
Collider. In theory, in order to improve on the accuracy of the QCD sum rule analysis for triply heavy 
baryons, one certainly needs to take into account the QCD 0{as) corrections to the sum rules in the 
further work. Additionally, it may be needed to carry out a comprehensive study on triply heavy baryon 
spectra from lattice QCD stimulations for the future. 
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